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TECHNICAL MEMORANDUM X-53257 


HEAT TRANSFER IN THE VICINITY OF A TWO-DIMENSIONAL PROTUBERANCE 

By 

Ed Murphy 

George C. Marshall Space Flight Center 
Huntsville, Alabama 

ABSTRACT 

Empirical equations are presented which will estimate the ratio of 
protuberance-to-flat plate heat transfer coefficients for supersonic 
flow in the vicinity of two-dimensional protuberances in a turbulent, 
supersonic boundary layer. The equations were derived from experimental 
heat transfer data taken around several two-dimensional protuberances, 
i.e., a 1 X 2-inch stringer, a 2 x 4- inch stringer with and without a 
30-degree wedge, and a one-quarter cylinder forebody. Important param- 
eters in the equations are boundary layer thickness, protuberance geometry, 
Reynolds number, and Mach number. Using the equations, the average devf- 
etlon between estimated and experimental ratios is from -9 
cent. 
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DEFINITION OF SYMBOLS 


Symbols 


h 

L 

L* 

M 

y 

5 

e 


Definition 

local convective heat transfer coefficient with 
protuberances 

local smooth flat plate convective heat transfer 
coefficient 

length from forward edge of protuberance 

length downstream of protuberance trailing edge 

Mach nuir±>er 

height of protuberance 

boundary layer thickness 

defined in Figure 1 
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HEAT TRANSFER IN THE VICINITY OF A TWO-DIMENSIONAL PROTUBERANCE 

SUMMARY 


Three empirical equations are presented for estimating the turbulent 
flat plate heat transfer coefficients in the vicinity of a two-dimensional 
surface protuberance under supersonic flow conditions. Each equation is 
applicable to one specific area in reference to the surface protuberance. 
The locations are the area forward of the protuberance, the separated 
flow area immediately behind the protuberance, and the area of reattached 
flow more than two protuberance heights behind the protuberance. The 
equations are limited to the vicinity of two-dimensional protuberances 
which have flat, cylindrically blunted, or wedge-type forward faces and 
a 90° flat rear face. The equations were developed from the experimental 
data of Burbank [3] using a dimensional analysis approach and reduced to 
their final form by trial and error. The increase in heat transfer coef- 
ficient is seen to be a function of boundary layer thickness, Mach number, 
protuberance geometry, and Reynolds number. A comparison between Burbank’s 
experimental data for several protuberance shapes with predicted values 
of h/ho from equations (1), (2), and (3) showed an average deviation of 
from -9 to +15 percent. 


I. INTRODUCTION 


In turbulent, supersonic flow, the local heat transfer coefficient 
for a vehicle skin with surface irregularities or protuberances can* vary 
widely from the local heat transfer coefficient for a skin which is 
unobstructed. For design purposes, it would be highly desirable to be 
able to calculate these "altered" heat transfer coefficients knowing 
the protuberance geometry and the "unaltered" air flow parameters such as 
local Mach number, Reynolds number, and boundary layer thickness, A 
recent literature survey showed that an analytical method for calculating 
local, supersonic heat transfer coefficients around protuberances was not 
in existence. However, Bertram [1] showed that, for hypersonic flow, 
empirical equations could be formulated to predict heat transfer coef- 
ficients around protuberances. Consequently, his method of formulation 
was extended to the supersonic regime. The heat transfer data around 
two-dimensional protuberances presented in Burbank’s report [3] provided 
excellent material for deriving the des ired ^empirical equations. 



II. TEST FACILITY AND MODELS 


The experimental data given in Burbank's report were obtained from 
P. S. Yip's protuberance heat transfer test [3] run in Langley Research 
Center's Unitary Plan Wind Tunnel. Test free stream Mach numbers were 
2.65, 3.51 and 4.44. Reynolds number per foot varied from 1.3 x 10® to 
4.7 X 10®. Three turbulent boundary layer thicknesses were achieved by 
mounting the test flat plate model at different locations within the 
tunnel test section. For thick boundary layers of six inches, the test 
plate was mounted on the tunnel side wall. For thinner boundary layer 
a splitter plate was installed in the test section. The instrumented 
surface was attached to the splitter plate, and the protuberances were 
mounted at different downstream locations giving boundary layer thick- 
nesses over the model of 0.7 inches and 1.5 inches. Boundary layer trips 
were used in the splitter plate configuration to assure turbulent flow. 

The two-dimensional protuberances tested are sketched in Figure 1. 
The basic models were a 1 x 2- inch stringer and a 2 x 4- inch stringer. 
Either a 30® wedge or a quarter-round cylinder could be added to the 
forward face of the 2 x 4-inch stringer in order to study the effects 
of leading edge geometry. A rather wide test plate and protuberances 
plus a centerline location of the heat transfer instrumentation were 
used in hopes of attaining and measuring two-dimensional flow. However, 
oil flow pictures taken at the end of the test hinted that true two- 
dimensional flow may not have been achieved. Therefore, the data used 
in deriving the two-dimensional protuberance equations may have some 
inherent inaccuracies caused by traces of three-dimensional flow charac- 
teristics in the flow field. 


III. DISCUSSION OF RESULTS 


As explained in the Appendix, dimensional analysis and a trial and 
error scheme to account for variable protuberance forward face configura- 
tions resulted in the following three equations for h/ho. Upstream of 
the protuberance, the heat transfer coefficients can be estimated by 
the equation 


h/ho = 60.0(1.10 


.10 5/y) 




( 1 ) 


if L/y S 4.0. The cot 0 term in the denominator provides a correction 


2 



for two-dimensional protuberances with wedge- type faces . Local Reynolds 
number is calculated from the free stream properties and the wetted 
length to the point being analyzed. The distance, L, is measured from 
the forward base of the protuberance to the point in question. Heat 
transfer coefficients behind the protuberance are approximated by two 
equations, the choice of which depends on the distance, L’ , downstream 
from the trailing surface. Immediately behind the protuberance where 
the flow is separated, the applicable equation is 


h/ho = I (L'/y)2 


1.0 + .021 



( 2 ) 


provided L*/y ^ 2, It was found that flow reattachment occurred at 
L*/y approximately equal to two regardless of the forebody shape. The 
following equation is used to estimate h/ho in the reattached area, 

L’ /y ^ 2.0. 


h/ho = 1.0 + .025 Re(y/L')2 M^(S/y)^/®. 


( 3 ) 


Notice that wake heating for two»^imensional protuberances is in no way 
affected by forebody angle. 

Figures 2 through 4 show the correlation of the function 





(y/s)^ 


from equation (1) to the ratio of protuberance- to-flat-plate heat trans- 
fer coefficients. Accuracy of the correlation is good for Mach numbers 
2.65 and 3.51. At Mach 4.44 the correlation is not nearly so good, 
mainly due to data scatter. Figure 5 shows the correlation of the 
function 


[Re (y/L/)^ M^(S/y)] 


of equation (3) with experimental h/h^ data. Again, data at Mach 4.44 
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show more deviation, The correlation of equation (2) with test data Is 
best seen In Figures 6 through 17. 

Figures 6 through 17 show profiles of the protuberance tested and 
the correlation of test data with all three equations. Reynolds number 
effects were sp small In the separated wake of each protuberance that 
only one calculated line was used In each figure. Accuracy of the cal- 
culated approximations is not only best for the lower Mach numbers, but 
also best for the more blunted protuberances. Equation (1) should be 
uaed wlthr caution for forebody, angles less rthan, 30 degrees, since accuracy 
seems to decrease with decreasing 6. Approximating a 45»degree forebody 
angle for the quarter-round cylinder forebody configuration (Figures IS 
through 17) yields reasonably accurate results. The most critical param- 
eters In these equations are Mach number, y/6, and y/L, Reynolds number 
effects are almost negligible and would be well within the accepted scatter 
of the original data, l.e., ±25 percent. 

Since equations (1), (2), and (3) were developed from data In the 
Mach number range 2.5 to 4.5, It would be wise not to extend their use 
to higher Mach numbers without first comparing them with experimental 
data. When used In this Mach number range, and for sufficiently blunt 
bodies, an accuracy of -9 to -flS percent can be expected In the predicted 
values of h/h^, 
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30® WEDGE 

FIGURE I. SKETCH OF PROTUBERANCE MODELS 
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FIGURE 2. CORRELATION OF EXPERIMENTAL DATA 

M = 2.65 
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FIGURE 3. CORRELATION OF EXPERIMENTAL DATA 

M = 3.51 
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FIGURE 4. CORRELATION OF EXPERIMENTAL DATA 

M = 4.44 








FIGURE 5. CORRELATION OF EXPERIMENTAL DATA 

EQUATION 3 PARAMETER 





Configuration- I x 2” Stringer 
M = 2.65 8 = 0.7 inches 
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FIGURE 6. COMPARISON OF EXPERIMENTAL TO EMPIRICAL PROTUBERANCE FACTORS 



Configuration- 1x2" Stringer 
M = 2.65 8 = 1.5 inches 
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FIGURE 7. COMmRISON OF EXPERIMENTAL TO EMPIRICAL PROTUBERANCE FACTORS 



Configuration* I x 2" Stringer 
M = 2.65 8 = 6.0 inches 

Re/ft X 10® Test Data 
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FIGURE 8. COMPARISON OF EXPERIMENTAL TO EMPIRICAL PROTUBERANCE FACTORS 
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FIGURE 9. COMPARISON OF EXPERIMENTAL TO EMPIRICAL PROTUBERANCE FACTORS 







Configuration* 2x4" + 30® Wedge 
M = 2.65 8 = 60 inches 
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FIGURE 12. COMPARISON OF EXPERIMENTAL TO EMPIRICAL PROTUBERANCE FACTORS 



Configuration- 2x4+ 30® Wedge 
M = 3.51 8 = 6.0 inches 
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FIGURE 13. COMPARISON OF EXPERIMENTAL TO EMPIRICAL PROTUBERANCE FACTORS 





Configuration* 2x4''+ 1/4 Round 
M s 2.65 8 * 6.0 inches 

Re/ft X 10® Test Data 
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FIGURE 15. COMPARISON OF EXPERIMENTAL TO EMPIRICAL PROTUBERANCE FACTORS 
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FIGURE 16. COMPARISON OF EXPERIMENTAL TO EMPIRICAL PROTUBERANCE FACTORS 



Configuration : 2 x 4" + 1/4 Round 
M s 4.44 8 * 6.0 inches 



21 


FIGURE 17. COMPARISON OF EXPERIMENTAL TO EMPIRICAL PROTUBERANCE FACTORS 
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FIGURE 18. DETERMINATION OF COEFFICIENT /9 FROM EQUATION 3 a 


22 








APPENDIX 

Method of Approach 


A study of the data from Reference 3 showed that the dimensionless 
protuberance factor, h/ho, depends on four independent variables: 


h/hQ = f(Re, y/6, L/y). 


(la) 


Assuming a product solution of the above parameters, an equation of the 
following form would result: 


h/ho = 0 Re^M^ (y/6)^ (Wy)^ (2a) 

The exponents a, 7, and | as well as the constant 0, can be detemined 
from the experimental data. 

The log of equation (2a) is 


ln(h/ho ~ In 0 + 0! In Re + p In + 7 ln(y/6) + | In (L/y). (3a) 


Taking the partial derivative of equation (3a) with respect to each of 
the variables, holding all the other parameters constant, allowed each 
exponent to be evaluated. For example 


Q In ^) 

o__ 

(b In M^) 


(In Re, In y/6. In L/y) 


constant 


= P 


where 0 can be easily determined by extracting the slope of the curve 
obtained by plotting h/ho versus on log-log paper as shown in Fig- 
ure 18. In this case, p is 3/4. The same operation, performed on all 
other variables, resulted in the following equation: 


h/h^ = 0 


m3 


-.1/4 


Re (L/y) (y/6)2 


(4a) 
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Equation (4a) was formulated using flat-faced protuberance data 
having 6/y ratios approximately equal to one. Data were available for 
various &/y ratios and also for protuberances with other forebody shapes 
such as a 30-degree wedge and a 1/4-round cylinder. 

By adding the factor y/2 cot 0 to L, equation (4a) would approximate 
the experimental data for the 30-degree wedge and quarter-round cylinder 
forebody cases. The y/2 cot 0 factor merely moved the point from which 
L was measured back one-half the distance to the original flange origin. 

The angle 0 for the quarter-round cylinder was considered to be 45 degrees. 

It was desired to extend the usefulness of the equation to include 
cases where the protuberance height was either less than, equal to, or 
greater than the boundary layer thickness. This was done by incorporating 
by trial and error a 6/y term into the equation. With the evaluation of 
the constant 0 from the test data, the final form of the equation is 


h/h = 60(1.10 - 
o 


.10 6/y) 



(y/6)2 J 


(5a) 


The equations for the other two regimes, the separated and reattached 
areas, were obtained in the same manner as equation (5a), with a similar 
correction being made for varying 6/y ratio. 
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